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We reconsider fluid dynamics for a self-propulsive swimmer in Stokes flow. With an exact 
definition of deformation of a swimmer, a proof is given to Purcell's scallop theorem including 
the body rotation. The breakdown of the theorem due to a finite Stokes number is discussed 
by using a perturbation expansion method and it is found that the breakdown generally occurs 
at the first order of the Stokes number. In addition, employing the Purcell's "scallop" model, 
we show that the theorem holds up to a higher order if the strokes of the swimmer has some 
^ , symmetry. 
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1 Introduction 

Fluid dynamics of locomotion of microorganisms such as bacteria and planktons has been 
studied for more than half a century ([E],[22]) and is still a hot topic in physics, mathematics 
and biology [H]. As a milestone to discuss the locomotion of such micros wimmers, there exists 
a well-known theorem called Purcell's scallop theorem [IS], which asserts that a microorganism 
with a reciprocal stroke in Stokes fluid cannot travel at all in one period of its motion. The 
proof of the theorem given in Purcell's famous lecture [19] was only schematic. Although 
the theorem were repeatedly discussed by many researchers ([3],[l],[7j,[l2],[20j) , the definition 
of the deformation of the swimmer appears not to have been paid much attention. Shapere 
and Wilczek [20J first established a theoretical formalism for a swimmer in viscous fluid in 
terms of gauge structure, and gave a proof for the scallop theorem. However, their theoretical 
framework based on gauge fields was conceptual and did not provide operational formulae for 
the locomotion of the swimmer. Yariv |23j improved their framework in order that we can 
calculate kinematic properties like velocity of the swimmer from its surface deformation from 
the fluid dynamical point of view. However, the scallop theorem was not completely proved in 
his paper, because according to his definition of the deformation of the body the rotation and 
the surface deformation of the swimmer are not uniquely identified as was pointed out by Yariv 
[23J himself. Also Childress and Dudley [5] states 

"As far as we know there has been no rigorous proof of this theorem based upon 

the mechanics of a Navier-Stokes fluid and free-swimming body" 

S_i 

on the present status of the scallop theorem. 

In the scallop theorem, on the other hand, the inertia of the fluid and the body are totally 

neglected. The scallop theorem and its breakdown have been recently looked back again [16J, 

and Childress & Dudley [5] discussed the possibility of sudden breakdown of the theorem at a 

nonzero critical Reynolds number. The breakdown of the theorem may also arise from multiple 

degrees of freedom of deformation and from fluid properties such as rheology [T3] and inertia 

[5]. [13]. The breakdown due to a finite mass of the swimmer, i.e. a finite Stokes number, 

was first discussed by Gonzalez- Rodriguez and Lauga [ID) . They provided general differential 

equations that govern locomotion of such a dense swimmer, and suggested that the scallop 

theorem does not hold at an arbitrary nonzero Stokes number. However, the definition of the 

deformation of a swimmer is not uniquely defined there, either. 
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Here, in this paper, we make clear the definition of the deformation of a swimmer, and give 
a rigorous proof for the scallop theorem including rotation together with translation. We then 
provide a perturbational argument on the breakdown of the scallop theorem by a finite inertial 
effect of a swimmer based on the equations of motion of both the fluid and the swimmer. 
Also with the Purcell's scallop model, we consider a higher order breakdown of the scallop 
theorem under some assumptions on a symmetry of the stroke of the swimmer. For discussion 
of these problems we develop a theoretical framework to describe the motion of microorganisms, 
introducing a virtual swimmer, which undergoes the same surface deformation as the real 
swimmer but without ambient fluid, to define rigorously the rotation and the deformation 
of the real swimmer. We call the coordinates attached to the virtual swimmer the vacuum 
coordinates, and the coordinates to the real swimmer the body coordinates. The transformation 
from the former to the latter defines the rotation (i.e. the gauge) of the swimmer. 

Some remarks should be made on physical situations of a finite Stokes number. As far as we 
consider typical microswimmers such as bacteria and mammalian sperms, the Reynolds number 
is small enough for the Stokes equation to be available. In reality, bodies of the microswimmers 
are usually a little heavier than the surrounding fluid. Among these microorganisms, relatively 
larger members like Volvox and Paramecium may have nonnegligible Stokes numbers, Rs ~ 1, 
and smaller Reynolds numbers, Re ~ 10~ 2 and Re ~ 10 _1 for each sample organism. A tiny 
bug in air may be another example for this situation because the averaged density of such a 
bug is much larger than that of the air. 

Here, in this paper, we consider the locomotion of the microorganisms which have a small 
but finite Stokes number in the fluid governed by the steady Stokes equation. Of course gravity 
effects on these swimmers become significant at high Stokes numbers, but in this paper we do 
not pay much attention to the gravity effects to keep the model as simple as possible. This paper 
consists of 5 sections. Section 1 is the introduction, and in section 2 we discuss a theoretical 
framework for a swimmer immersed in Stokes fluid and give an equation that governs the 
locomotion of the swimmer. In section 3 we restate the scallop theorem and give a complete 
proof to the theorem. In section 4, a perturbational discussion of the breakdown of the scallop 
theorem due to a finite Stokes number is made for a swimmer without rotation. In section 5 
we consider the breakdown of the scallop theorem under the assumption of a symmetry of the 
stroke using the Purcell's scallop model [TJ. Summary and conclusion are given in section 6. 

2 Formulation 

In this section we set up formulae for velocity and angular velocity of a self-propulsive swimmer 
in fluid. To discuss the motion of the swimmer, we first define a virtual swimmer, which deforms 
its body in exactly the same way as the real swimmer except that the virtual swimmer has 
no surrounding fluid, and thus experiences no external forces. The virtual swimmer therefore 
conserves the total momentum and the total angular momentum both of which we assume 
to be zero. We further assume that at the initial time the virtual swimmer exactly coincides 
with the real swimmer, and therefore both their centers of mass locate at the same position, 
and their orientations are the same. We attach an inertial coordinates, which we call the 
vacuum coordinates, to this virtual swimmer with its origin located at the center of mass of 
the virtual swimmer. The virtual swimmer has exactly the same shape as the real swimmer 
at the same time, and thus the latter is obtained by an affine transformation from the former. 
We also define the body coordinates attached to the real swimmer as the coordinates obtained 
by the same affine transformation from the vacuum coordinates. Thus the origin of the body 
coordinates is located at the center of mass of the real swimmer. We denote the orthonormal 



basis of the vacuum coordinates by e»(i = 1, 2, 3), which is independent of time (Fig. [Q. The 
motion of the virtual swimmer is described in Lagrangian coordinates where a position of a 
Lagrangian particle of the virtual swimmer, f(a,t) = Ylifi( a it)ei, is regarded as a function 
of the Lagrangian coordinates a = (ai, a 2 , 0,3) and time t with f(a, 0) = a . 

We now define the surface deformation velocity u' of the virtual swimmer as u' = ^2 i -^-(o,, t)ei 
where a is assumed to be on the surface of the swimmer. The real swimmer, however, not only 
deforms but also translates and rotates under the action of the external force from the surround- 
ing fluid. The translation velocity U is defined as U = dX/dt where X(t) is the center of mass 
of the real swimmer. The position of the Lagrangian particle of the real swimmer f(a, t) with 
respect to the center of mass of the real swimmer is then obtained as f(a,t) = R(t)f(a,t), 
where R(t) is a rotation matrix in SO (3) and -R(O) = 1, because the real and virtual swimmers 
have exactly the same shape. The unit vectors of the body coordinates e\(i = 1,2,3) is then 
obtained as Bi(t) = R(t)ei. 

The total surface velocity u of the real swimmer is 



u 






u + n(t) x f + u', 



(2.1) 



where we have written / as / = Xl«/* e *' an d used the fact that fi(a,t) = fi(a,t) and e-i = 
fixef The last term u' means the surface deformation velocity of the real swimmer defined 
as it' = Y^i ~~it&i which is equal to Ru' . The rotational angular velocity vector fl(t) is defined 
together with two skew symmetric matrices A(t), B(t) as 



dR(t) 

dt 

MM 



B(t)R(t) = R(t)A(t) = fix R(t), 



(2.2) 
(2.3) 



where the Einstein convention for repeated indices is employed. The angular velocity vector 
O(t) can also be written as fl(t) = (1/2) £\ e, x gj. 
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Fig. 1: The vacuum coordinates {ej} and the body coordinates {ei}. 



The governing equation of the fluid surrounding the real swimmer is the incompressible 
Navier-Stokes equation. The equation for fluid velocity v = v(x,t) is written in a non- 



dimensional form as follows: 

dv 
V-a = R U] —+Re(v-V)v (2.4) 

at 

V-v = (2.5) 

and the stress tensor er is given by 

tr = -pl+ (Vv + (Vv) T ) . (2.6) 

The non-dimensional parameters of the equation, Re and R^, are the Reynolds number Re = 
pVL/p and the oscillatory Reynolds number R w = pL 2 u/p, the latter of which corresponds to 
the product of the Reynolds number and the Strouhal number St = Luj/V , where L, V and 
cu are respectively the characteristic scales of length, velocity and frequency (d/dt) of the fluid 
motion, with the density p and the viscosity p of the fluid being assumed constant in this paper. 
The boundary condition for v at the body surface of the real swimmer is that v coincides with 
the total surface velocity u of the real swimmer, while v — y at infinity. 

From Newton's equation of motion for the real swimmer, we obtain the following equation 
in the non-dimensional form, 

*s(m) = ©- (27) 

where J is the inertial momentum tensor of the real swimmer. The non-dimensional number 
Rs, the Stokes number, is represented as R$ = PmL 2 oj / p where pu is the mean density of the 
swimmer. The force F and the torque T acting on the real swimmer also non-dimensionalized 
like F* = pLVF using the Stokes law of resistance, where the asterisk denotes the dimensional 
quantity. This non-dimensionalization is effective in a situation without other external forces 
such as gravity and electromagnetic forces. 

Hereafter, we assume that the non-dimensional parameters ReC R^CRs satisfy the inequal- 
ity, Re, Ru <C Rs <C 1. In this section we consider the motion of the swimmer in the fluid with 
Re = R w = 0. We will derive equations for the velocity U and the angular velocity ft. When 
the fluid obeys the steady Stokes equation, Lorentz' reciprocal theorem [18] gives 



dSn ■ cru = / dSn ■ &u, (2.8) 

s Js 

where S denotes the surface of the real swimmer, and the symbol "hat" indicates quantities 
of another solution of the Stokes equations with the same boundary shape S and the vector n 
denotes the unit normal vector to the surface. For the solution u we take a solution satisfying 
the boundary condition, 

u = U + Qx f (2.9) 

at the surface S of the swimmer [21] . With linearity of the Stokes equation, the stress tensor 
is given by 

<r = t T U + t R Q, (2.10) 

where S^ and S/j are third rank tensors [UJ which depend upon the direction of the body 
coordinates e^. Substituting equations (12.91) and (12.101) into (12. 8p we have 

fw ^ - f *s (t u :»v (*) , cm) 



TJ \nj J s \(n-'E R ) T uJ \n 



where the superscript T denotes the transpose of a matrix. Here n • Sy and n ■ S^ denote 
the second rank tensors nJlTijk and niEmjk respectively. Using the 6x6 symmetric resistive 
matrix Q| 

'K T , K£\ _ f, Q f n-XJr, n-T, R 



K =u;; £j=y/H/x(™.^), /*<».£,)) < 2i2 > 



and arbitrariness of £/ and fi, we obtain 



^ tffl+ftsit&Z)- (»3) 



With Newton's equation of motion (J2.7)) . we obtain the desired equation 

d i U \ I Ftr\ I Fd 



»**\m) = Ktt R ) + V°) (214) 

S ^ * © (215) 

F D \ f JO f(n-i: T ) T u' 



r » , dS( K - 1 T r, . (2.16) 

The force .FVi? and the torque T Tfi are the external force and the torque from the fluid, 
arising from the translation and rotation of the real swimmer, while the force Fp and the 
torque Tp arise from the surface deformation of the swimmer. We remark that the right hand 
side of equation (12. 14ft represents the total force and torque exerting on the real swimmer, and 
therefore when we take into account the gravity and the buoyancy effects, we only need to add 
the gravity and the buoyancy forces to F and their torques to T, respectively. We should note 
that both the third rank tensors, Sy and S#, and the resistive matrix K of the real swimmer 
depend only on the surface shape of the real swimmer /(£). 

As we are considering the scallop theorem which is concerned with the motion of the swim- 
mer due to its surface deformation defined by using the virtual swimmer, it is convenient to 
describe the problem in terms of quantities of the virtual swimmer in the vacuum coordinates. 
We then have the following equations; 

/tt\ r //„ -^ \T„./\ 

(2.17) 



d ( U \ 

Rs dl {RiVR-'n) - RKR 


*©W.«(f::K 


where the 6x6 matrix R is defined as 




R = 


(-«)• 



(2.18) 

1 The relation of off-diagonal components of K is as below using Lorentz' reciprocal theorem. If we take two 
solutions with the boundary condition u = U and u = fix f respectively at the surface S 1 , we obtain er = E^C/ 
and <r = E^Q . Substitution of these into equation (|2.8[) gives 



/ dSn ■ (E T U)n x / = / dSn ■ (E R il)U. 
Js Js 

After some manipulations, we get 

/ dSn ■ f x (n • Er)*/ = / dStl ■ (n • E fl ) T £/, 
Js Js 

which implies that the transpose of an off-diagonal component is another off-diagonal component. 



and I v is the inertial momentum tensor of the virtual swimmer. When the deformation velocity 
of the virtual swimmer is given, this equation together with ( 12. 2 j) determines the rotation matrix 
R(t) and the translational velocity U(t). 

3 The proof of the scallop theorem 

In this section we give a proof of the scallop theorem using equations (I2.2p and ( 12.1 7p in the 
case of the vanishing Stokes number, R$ = 0. 

We consider the case where the shape of the virtual swimmer deforms in a reciprocal manner, 
i.e. the shape once deformed retraces back to the initial shape. The mathematical definition of 
the reciprocal motion is that for the surface deformation of the virtual swimmer, which starts 
at t — and ends at t = T, there exists a continuous function g(t) such that f(t) = Q(t)f(g(t)) 
and g(0) = g{T) = where Q(t) € 50(3) is a three dimensional rotation matrix, allowing the 
possibility that the swimmer takes different directions at time t and g{t). We assume that g(t) 
is smooth except at a finite number of points, and then intervals of integration over the time 
t we have in the following in this paper should be divided into those in which g(t) remains 
smooth. In this case, we can prove that Q(t) = 1 as below. 

The total angular momentum of the virtual swimmer always vanishes, and therefore 

f df f dQf(a,g(t)) 
= / Prnf{a,t) x —(a,t)da = / p m Qf(a,g(t)) x — da. (3.1) 

Noting that 

= Oft x f(a, t) + -Q-(a, g(t)), (3.2) 

where the angular velocity vector O^ is defined as 

f^x Q , n^-e^fo-)^, (3.3) 

and p m = Pm{. a ) is the density of the swimmer, we find the contribution from the second term 
of equation (13.21) to (13. ip vanishes as 

p m Q(t)f(a, g(t))x ^Q%(a, g(t)) da = ^Q(t) J p m f(a, g(t))x ^(a, g(t)) da = (3.4) 

due to the vanishing initial angular momentum conserved. Then the equation (I3.2p gives 

J p m f(a, t) x (n Q (t) x f(a, t))da = I v (t)n Q (t) = 0. (3.5) 

For a general 3-dimensional swimmer, I v is not degenerated and we have Q®(t) = 0. Then 
equation (13. 3p with Q(0) = 1 gives Q(t) = 1. 

The scallop theorem asserts that the position and the direction of the real swimmer at the 
final time t = T coincide with the initial position and direction if the motion of the swimmer 
is reciprocal, the surrounding fluid obeys the steady Stokes equation, and the Stokes number 
of the swimmer vanishes. In this case, the equation (I2.17P is reduced to 



where the 6x6 matrix M, the mobility matrix, is the inverse of the resistive matrix K. 
Denoting if = g(t) for short, the equation (13. 6p is reduced to 

u (t)\ . D^u/im (hit) dS ( n ■ ^T) T (t) ■ m/dtW 



^r-^H/^i-w-i^ 1 (3 - 7) 

R(tM(t>) I V) dS ^ n ■ ^^ • P/^M ^ f3 g , 

) _ l(f/i /J7(t')\^ / 



where we have used the fact that M(t) = M(t'), 5(t) = S(tf), S T (t) = £ T (t') and E fl (t) = 
E#(£'), which holds because f(a,t) = f(a,t'). Especially from this equation, we obtain 

Hi' 

R- l {t)Q{t) = R- l {t')Q{t')^-. (3.10) 

Using equations (12. 2 j) and (12. 3 p we obtain 

^j(i) = -Rj/ B ik R k j = —€ikpRnRkjQp, (3-11) 

and from the definition of the determinant of the matrix i? we have 

which leads to 

^Ikr-^-li^kj ^ijq-^rq ^ijqyrL )q r yO.i-O) 

by the use of \R\ — 1 and the multiplication of R rq to equation (I3.12p . Substituting equation 
(I3.13P into equation ( 13. lip , we have 

A tJ (t) = -e ijk (R- 1 (t)n(t)) k , (3.14) 

which together with (13.101) gives 

A(t) = A(t')^. (3.15) 



The rotation matrix R(t) satisfies 

dR{t) 



R(t)A(t)., (3.16) 



dt 
and R(g(t)) is also the solution of (I3.16p . because 

^» = ^f faW) = !fW«)W»<«» = WW) 0.17) 

and H(#(0)) = Jl(0) = 1. Therefore #(#(£)) = i?(t), and thus R(g(T)) = R(0) = 10. 
Using this relation in equation (13. 9p . we have U(i) = U(t')dt'/dt which implies 

d _x^ = dj^dx = *m = 

dt dt dt ViA ;; dt vyw; v ; v ; 

Together with X{g{0)) = X{0), this means that X{g{t)) = X{t) and thus X(T) = X(g(T)) = 
X(0) = 0, which completes the proof of the scallop theorem. 



2 This can be proved also by noticing that 

R{T) = R(0)Te& A{t)dt = R(0)Te^ A ^ dt ' = 1. 
where L symbolize the round integration from t' = to t' — 0, and T means anti time-ordering operator. 



4 Breakdown of the scallop theorem due to a finite mass of the swimmer 

Let us discuss the breakdown of the scallop theorem by a nonzero Stokes number. Assuming 
that the swimmer moves in the fluid with Re = R^ = 0, and the Stokes number is a small but 
nonzero constant, we employ -R^-expansion as R, X, Q and U: 

R = R (l + R s Ri + ---) (4.1) 

X = X (0) + R S X {1) + R 2 S X {2) + ■■■ (4.2) 

n = ra (0) + R s n w + R 2 s n {2) + ■■■ (4.3) 

U = C/ (0) + RsZJM + i?2[/(2) + . . . _ (4.4) 

At the first order of the Stokes number, substitution of these expansions into (12.21) and (12.171) 
gives 

Q« xR = R —. 4.5 

at 

and 

5 (jwiu' • n<»>) = R ° (RlK - KRl)R » M + R ° KR ° (n«) • < 4 ' 6 > 

where the matrices Ro and Ri are 6x6 matrices defined as 

Here we focus our attention on a swimmer without rotational motion (Q(t) = 0). Then we 
can assume that the motion of the swimmer is in the direction of e 1; reducing the equation 
(I2TT7D into 

dU(t) 

e-±l=-K(t)U + F(t), (4.8) 

where we denote e = Rs, U(t) = U(t) ■ e\, K(t) = —(K T (t)ei) ■ e 1; and F(t) = F D (t) ■ e l7 and 
K(t) is positive definite. Then the equation (14 .8p is exactly solved as 

U(t) = l e -$SSKV)# f d t'F(t>)e^° K ^ dt " + U(0)e'^o K ^ dt '. (4.9) 

e Jo 

Integrating by parts repeatedljo, we represent the equation (14. 9 j) as the form of the asymp- 
totic power series of the Stokes number e as U(t) ~ Y^^=o £ n U( n '(t), where we use an assumption 
that K(t) and F(t) are smooth functions, and eliminate exponentially small terms from the 
asymptotic expansion. Then the first few terms of the expansion are 



U^(t) = -j t {F(t)K-\t))K-\t) (4.11) 

UW ® = T+ \ i{mK-\t)) K-\t)) K-\t). (4.12) 



^Denoting E(i) — e? J« Kiyt ^ dt , we have E{t) = eK 1 {t)E(t) and the integration becomes 
F{t')E{t')dt' = e [E{t')F{t')K-\t%~e f A (F{t')K-\t')) E(t')dt> 
and similarly. 



o Jo dt 



We now consider the breakdown of the scallop theorem. Assume a swimmer in a reciprocal 



motion. Since the propulsive force F satisfies F(t) = F(t')% and K(t) = K(t'), where t' = g(t), 
the zeroth order of equation (14. 8 p gives 

X (°)(T)= [ F(t)K-\t)dt= [ Ftf)K-\lf)dlf = 0, (4.13) 

Jo Jo 

which corresponds to the scallop theorenoD According to equation (14. lip , the first-order dis- 
placement in one period of the reciprocal motion is represented as 

r T r T du {0) (t) 

X (1) = / dtU{t) = - dt — -±±K-\t). (4.14) 

Jo Jo dt 

In the case of the reciprocal motion as g(t) = t on [0, T r ] and g(t) = (T — t)T r /(T — T r ) on 
[T r ,T fl, equation (T4TT4"P leads to 

m T r Tr , n s dK- l (t) 

A (1) (T) = 7F ^ F - / dtU {0 \t)- 



T — T r J dt 

+ — ^- ([/(°)(+0)^ 1 (+0) - U<®{T r - 0)K-\T r - 0)) . (4.15) 

This does not necessarily vanish as we show in the next section That X^ does not generally 
vanish, means that the scallop theorem breaks down at the first order of the Stokes number. 
This results are consistent with the continuous breakdown of the scallop theorem suggested by 
Gonzalez- Rodriguez and Lauga [10]. If, however, the stroke of the swimmer has some symmetry, 
the first-order displacement in one period of the stroke can also vanish, and the theorem holds 
true up to a higher order of the Stokes number, as also discussed in the next section. 

5 Purcell's scallop model 

In this section we discuss the motion of the Purcell's "scallop" model [1] as an example of a 
swimmer with a finite Stokes number (Fig. [2]). This swimmer with one hinge to move was first 
introduced by Purcell to explain the scallop theorem [TS]. The "scallop" has two slender rods 
with the same length / and with the same circular cross section of radius r. From the symmetry 
of the shape, it does not rotate and so become an example of the swimmer without rotationD 
We calculate resistive and propulsive forces using Cox's slender-body theory [6], which describes 
approximately well when the aspect ratio of the body a = r/l is much smaller than 1. Let us 
take a coordinate si— I < s < Zjalong the body from the bottom edge to the top, with the point 
of s = at the hinge. Let X(s) denote the normal vector to the cross sectionD Cox's theory 
gives the force exerted by fluid on an infinitesimal part of the body between s and s + ds as 

dF(s)=( n-a-dl) ds = -((2- XX) -uds, (5.1) 

where ( = 27r/i(log (2/a))" 1 . Here the line integral is taken along the circumference L(s) of the 
cross section at the position s, and u(s) indicates the total surface velocity of the part of the 
body, which is equivalent to the fluid velocity at the surface. 



4 Here again we have used the symbolic notation L . 

5 The dg(t)/dt is discontinuous at t = T/2, and we divide the integral interval into [0, T r ] and [T r ,r] 




Fig. 2: The conceptual figure of the virtual swimmer of the Purcell's "scallop" D The coordinates 
ei, B2 are the vacuum coordinates with the origin at the center of mass of the virtual "scallop" . 



The scallop is configured in two dimensional e^-plane as shown in Fig. [2j The motion 
is then in one- dimensional direction. The motion is driven by temporal variation of the angle 
a; a(t) changes reciprocally as a — > a m — > a (Fig. [3]). The slender-body theory gives 
Il(s) dl{n ■ St) t and the resistive coefficient K as 



f dl(n--E T ) T (s) 

JL{s) 



K 



< 



2(1(2 



2 — cos 2 a 
-sgn(s) cos a sin a 



-sgn(sj cos asm a 
2 „;„2„ I 



sin 2 a 



cos a) 



(5.2) 
(5.3) 



in the vacuum coordinates. When we write the position and the velocity of the hinge as 
Xnif) = —0.5/ cos a and £/#(£) = +0.51 sin a (da /dt), it is found that the surface deformation 
velocity of the virtual swimmer is u'(s) = (Ujj, 0) + s(d\/dt), where a = a(t) depends on time 
t and hereafter we assume some smoothness of a(i) so that d(0) = a(0) = 0. The propulsive 
force F is then derived from (15.21) as 



F = 2(l 2 sin a cos 2 a 



da 
~dt' 



Substitution of (15. 3p and (15. 4p into equation (14. 8 p gives 



,dU 
~dt 



9 9 CtC* 

-(2 — cos a)U + cos a sin ex— I, 

(II 



(5.4) 



(5.5) 



where e' is defined as e' = a 2 log(2/a)e/2cj, and u is the frequency of the stroke of the swimmer. 
Henceforth we nondimensionalize quantities, using I and 1/u for the unit of length and time. 
Let us assume e'<l and consider the expansion of the equation (15.51) in terms of e'D At the 
zeroth order of e'Cwe get the velocity U^ as 



U^(t) 



cos a sin a 



cos^ a 



da 
~dt 



(5.6) 



and thus the displacement of the center of mass in one period becomes 



l(°>(t)= / U^(t')dt' 



a cos 2 a' sin a' 



oil 



2 — cos 2 a' 



■da'. 



(5.7) 
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Since X^ is a function of a, the realization of the scallop theorem is confirmed. At the first 
order of e', equation (14. lip gives the first order displacement as 



X«(T) 



UW[t)dt = - 



T 2 sin 2 a cos 3 a ( da 



(2 - cos 2 a) 2 \ dt 



dt. 



(5.8) 



Here we have used integration by parts and d(0) =0. 

We now assume a reciprocal motion in which the swimmer opens and then closes its "shell" 
during one period. If < a < a < a m < n/2, then the integrand of equation (15.81) is always 
positive, and the swimmer moves to —e.\ direction by the effect of the mass inertia of the 
swimmer. 




Fig. 3: Temporal variation of the shape of the swimmer. 



If < «o < at < a m < 77 and a r 



7[ 



«o, the period of motion of the swimmer can then 



be divided into 4 intervals, I ao ->-K/2i In/2->a m , Ia m ->-K/2i Iir/2^a , m eacn °f which the motion is 
time-reversal or mirror-symmetric to that in another interval. If the speed of the deformation 



in I. 



X, 



(i) 



7r/2-S>Q m 



is C2 times faster than that in I a ^. w / 2 , the first-order displacement in the first half, 



oo— >a m i 



becomes 






;i-c 2 )x 



(1) 

QO— >7r/2' 



(5.9) 



where X 



(i) 

ao—nr/2 



denotes the first-order displacement in I ao ^-x/2- Similarly, if the deformation 
speeds of in J am _ >7r /2 and Iir/2-^ao are by C3 and C4 times faster than that in I ao ->n/2 respectively, 
the net displacement in one period of the reciprocal motion is 



X«(T) = (l-c 2 -c 3 + c 4 )X, 



(i) 

dQ— >tt/2' 



(5.10) 



If 1 — C2 — C3 + C4 = 0, we obtain X^(T) = and the scallop theorem still holds up to the first 
order of Stokes number. 

In the second order of e', we have 



X( 2 )(T) 



1 



2 — cos 2 a 



'dU^ 
~dT 



dt 



T 4 sin 3 a cos 4 a 
(2 — cos 2 a) 5 




(5.11) 



using integration by parts and the assumption that the angle a = a(t) is a smooth enough 
function of time to satisfy a(0) = 0. With the same assumptions on the speed of the motion 
as in equation (15.101) . the second-order displacement is obtained as, 



X&>(T) 



[l + cl-cl- c\)X^ 



ao—yir/2' 



(5.12) 
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Thus the symmetric deformation of (c 2 , c 3 , c 4 ) = (1, 1, 1) gives the zero net displacement up to 
the second order: X^'(T) = X^ 2 '{T) = 0. On the other hand, the different deformation speeds 
as (c2, C3, C4) = (1, 2, 2) produces the net displacement X^{T) 7^ in one period of the motion 
while the first order displacement then vanishes: X^ l '{T) = 0. This result implies that if the 
"shell" opens slowly and closes quickly, the swimmer with the symmetric stroke shown in Fig. 
Ogives the net displacement in — e.\ direction at the second order of Stokes number. 

Before ending this section, let us consider the case of £ 3> 1. Using equation ( 14.91) . we 
expand the velocity of the swimmer in Taylor series of 1/e: 



t rt ft' 



U(t) = (l/e) dt'F(t') + (l/e 2 ) dt'F(t') K(t")dt" + 0{lje 6 ). (5.13) 

Jo Jo Jt 

Applying this equation to equation (15.51) . we obtain at the first order of 1/e' 

U {l \t)=j cos 2 asinada, (5.14) 

which is independent of the speed of deformation. In the case of < a^ < a < a rn < ir/2, X^ 
is positive quantity. Thus the "scallop" goes to +ei direction when its "shell" opens slowly 
and closes quickly as real scallops do. 

6 Summary and Conclusion 

We have established a framework to discuss a motion of a swimmer immersed in Stokes fluid, 
and given a rigorous proof for the scallop theorem. We have also discussed the breakdown of the 
theorem due to a nonzero mass of the swimmer, and shown that the degree of the breakdown 
depends on the symmetry of the stroke using the Purcell's scallop model. 

First of all, in order to define the deformation velocity we introduced the virtual swimmer 
which has the same shape as the real swimmer but has no ambient fluid. We then attached the 
the vacuum coordinates to the virtual swimmer and the body coordinates to the real swimmer. 
The position and orientation of the real swimmer is obtained by the affine transformation 
from the virtual swimmer. We derived the formulae which provide the velocity and the angular 
velocity of the swimmer when the surface deformation of the virtual swimmer is given in the case 
of vanishing Reynolds number Re and oscillatory Reynolds number R^. Using these formulae, 
we have proved the scallop theorem at a vanishing mass of the swimmer, or equivalently a zero 
Stokes number. 

Then we studied the breakdown of the scallop theorem, taking a finite Stokes number 
into consideration especially in the case of the swimmer without rotating motion. We showed 
the net displacement is generally at the first order of the Stokes number using an asymptotic 
expansion. We took Purcell's scallop model with a finite mass as an example and demonstrated 
the breakdown of the scallop theorem, but also showed that the theorem holds up to a higher 
order of the Stokes number for a swimmer deforming with a particular symmetry. 

Our argument is based on the assumption of Re, R w <C 1 and Re, R w <C Rs- In reality, 
however, the time derivative term of the fluid equation may become important for such a dense 
swimmer, where Re <C R^,Rs, and ife< 1 are satisfied. Gavze [9] gave formulae for the 
fluid dynamical force and torque on a rigid (non-deformable) body with an arbitrary shape 
under such conditions. However, to our knowledge, an explicit form of fluid dynamical force 
on self-propulsive deformable swimmer has not yet been found in non-stationary Stokes flow. 
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It is still an open question whether the scallop theorem holds in the case of a finite oscillatory 
Reynolds number. 

The gravity effect on a swimmer may also be important. In the case of Volvox, the mass 
density of the body is a little heavier than the surrounding water. According to Drescher et 
al. [8] , the density difference between the Volvox and water is approximately Ap ~ 2 x 10 _3 g/cm 3 . 
The characteristic velocity due to the gravity effect is estimated to be U g ~ 2 x 10 2 cm/sec by 
balancing the gravity effect with the resistive force on the body, 

lirA'Apg i 



6nfxAU g 

This velocity is comparable to the propulsive velocity of the swimmer. Burton et al. [2] has 
recently studied a neutrally buoyant Purcell's scallop model with separated positions of the 
centers of mass and of buoyancy. This kind of separation often occurs in microorganisms due 
to heterogeneity of mass distribution, and it may be of interest to apply our formulation to 
these organisms by including gravity and buoyancy. 

Before ending, we should remark the possibility of another choice of the gauge fixing to 
define the deformation velocity of a swimmer. In this paper we have introduced a virtual 
swimmer and its associated vacuum coordinates to define the deformation velocity. This choice 
appears natural but is not unique, and there may be another coordinate system, which is useful 
in considering other subjects beyond the scallop theorem. 
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